note from lecture: Multinomial regression coefficients indicate the change in log of relative risk ratio (conditional log odds) relative to base category. Can interpret direction, significance, and relative magnitude, but hard to evaluate magnitude.
A one-year increase in education is associated with a 0.15 unit increase in the conditional log odds of believing the bible is the inspired word of God, but not to be taken literally. In other words, holding gender and religion constant, people with more years of education are more likely to believe that the bible is the inspired word of God, rather than to believe it is the “Word of God.”
A one-year increase in education is associated with a 0.19 unit increase in the conditional log odds of believing the bible is an ancient book. In other words, holding gender and religion constant, people with more years of education are more likely to believe that the bible is an ancient book, rather than to believe it is the “Word of God.”
1B
Holding gender and education constant,
Those who reported having no religious identity are most likely to believe the Bible is an ancient book relative to believing it is the “Word of God.”
Those who reported having Jewish identity are most likely to believe the Bible is the “word of God” relative to believing it is the “inspired word.”
2A
I am using data from the Health and Retirement Survey, a nationally-representative survey of U.S. older adults. The sample consists of respondents to the survey’s subsection on Widowhood and Divorce, who experienced the death of a spouse between 2018 and 2022. This subsection consists of questions about the financial impacts of the death, such as changes in income, social assistance and work hours, changes in insurance coverage, death expenses, and so on.
Dependent variable:
Source of financing funerals: Did the deceased’s insurance and estate fully cover death expenses (a.k.a. the widow reports paying less than $300 beyond the covered amount), or did the widow have to rely on selling assets & withdrawing savings, relatives & friends, charity institutions, loans, or other sources?
Independent variables:
Widow’s demographic characteristics: Gender, race & ethnicity, foreign-born status
Widow’s employment status: currently employed or not (note: “non-employed” includes unemployed, retired, homemaker, disabled)
2B
library(pacman)p_load(tidyverse, broom, haven, skimr, janitor, marginaleffects, lmtest, modelsummary, flextable)# nnet package for multinomial logit# install.packages("nnet")library(nnet)# Turn off scientific notationoptions(scipen =100)# refresh environmentrm(list =ls()) # set working directory to project directory# setwd(here::here("FINAL_PROJECT_HRS_SPOUSAL_DEATH"))# read data in relation to working directorydf <-read.csv("../Input/HRS_widows_employ_tracker_cleaned.csv")# set thousand dollar units to aid interpretationdf$deathexpense_1k <- df$deathexpense_usd/1000# convert variable to factor and set base category df$deathexpense_sources <-as.factor(df$deathexpense_sources)df$deathexpense_sources <-relevel(df$deathexpense_sources, ref ="insurance_estate_full")base_category <-levels(df$deathexpense_sources)[1]print(base_category)
# weights: 48 (35 variable)
initial value 632.491093
iter 10 value 467.546076
iter 20 value 458.309198
iter 30 value 457.314866
iter 40 value 457.278160
final value 457.278087
converged
The base category here stands for cases where the widow had to sell assets or withdraw savings to cover death expenses.
model <-multinom(deathexpense_sources ~ foreign + female + employed + black_nh + hisp_allraces + non_bwh, data = df)
# weights: 48 (35 variable)
initial value 632.491093
iter 10 value 471.682568
iter 20 value 458.342032
iter 30 value 457.341989
iter 40 value 457.278218
final value 457.278087
converged
# weights: 48 (35 variable)
initial value 632.491093
iter 10 value 467.546076
iter 20 value 458.309198
iter 30 value 457.314866
iter 40 value 457.278160
final value 457.278087
converged
# Save AME output ame1 <-avg_slopes(model)# Make table with AME outputmodelsummary(ame1, shape = term : contrast ~ group,stars = T,notes =str_c("N = ", glance(model)$nobs),title ="Multinomial Logit")
tinytable_b36zz5olc9u7e5f8t3do
(1)
Multinomial Logit
assets_savings
charities
insurance_estate_full
loans
other
relatives_friends
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001
N = 353
black_nh 1 - 0
-0.092
0.023
-0.023
-0.015**
-0.074**
0.182**
(0.061)
(0.030)
(0.067)
(0.006)
(0.026)
(0.067)
employed 1 - 0
0.136*
-0.024
-0.120*
0.004
0.055
-0.051
(0.063)
(0.017)
(0.059)
(0.015)
(0.041)
(0.046)
female 1 - 0
0.090+
-0.005
-0.040
-0.015
-0.042
0.013
(0.051)
(0.020)
(0.055)
(0.020)
(0.036)
(0.045)
foreign 1 - 0
0.195*
0.094
-0.262***
0.003
0.038
-0.068
(0.090)
(0.072)
(0.057)
(0.013)
(0.047)
(0.050)
hisp_allraces 1 - 0
-0.345***
-0.035*
-0.036
0.077+
0.027
0.312***
(0.042)
(0.017)
(0.082)
(0.043)
(0.044)
(0.084)
non_bwh 1 - 0
-0.122
-0.033***
0.152
-0.012*
0.020
-0.005
(0.087)
(0.010)
(0.124)
(0.005)
(0.071)
(0.112)
Num.Obs.
353
R2
0.253
R2 Adj.
0.252
AIC
984.6
BIC
1119.9
RMSE
0.34
It seems that Hispanic identity has the most impact on the probability that someone relies on relatives and friend (31% higher probability than non-Hispanics) and on assets and savings (35% lower probability than non-Hispanics).
2E
plot_predictions(model,# the second argument must be "group"condition =c("foreign", "group"))
2F
Holding all other variables constant, to cover their spouse’s death expenses,
Foreign-born widows are significantly less likely than US-born widows to have their spouse’s death expenses fully covered by their spouse’s insurance or estate.
Hispanic widows are significantly less likely than non-Hispanic widows to relieve themselves of assets & savings and significantly more likely to rely on relatives and friends.